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ABSTRACT 

We report observations of dust continuum emission at 1.2 mm towards the star forming region NGC 
6334 made with the SEST SIMBA bolometer array. The observations cover an area of ~ 2 square 
degrees with approximately uniform noise. We detected 181 clumps spanning almost three orders of 
magnitude in mass (3M — 6 x 10 3 M Q ) and with sizes in the range 0.1-1.0 pc. We find that the clump 
mass function dN/dlog M is well fit with a power law of the mass with exponent —0.6 (or equivalently 
dN/dM oc M -1,6 ). The derived exponent is similar to those obtained from molecular line emission 
surveys and is significantly different from that of the stellar initial mass function. We investigated 
changes in the mass spectrum by changing the assumptions on the temperature distribution of the 
clumps and on the contribution of free- free emission to the 1.2 mm emission, and found little changes 
on the exponent. The Cumulative Mass Distribution Function is also analyzed giving consistent results 
in a mass range excluding the high-mass end where a power-law fit is no longer valid. The masses and 
sizes of the clumps observed in NGC 6334 indicate that they are not direct progenitors of stars and 
that the process of fragmentation determines the distribution of masses later on or occurs at smaller 
spatial scales. 

The spatial distribution of the clumps in NGC 6334 reveals clustering which is strikingly similar to 
that exhibited by young stars in other star forming regions. A power law fit to the surface density of 
companions gives £ oc #~ 062 . 

Subject headings: stars: formation — molecular clouds: individual(NGC 6334) — stars: initial mass 
function 



1. INTRODUCTION 

1.1. Massive Star Formation and Molecular Cloud 
Structure 

Half of the mass in the interstellar medium is in the 
form of molecular gas exhibiting a broad range of struc- 
tures, ranging from small isolated clouds, with masses of 
a few Mq and subparsec sizes, to Giant Molecular Clouds 
(GMCs), with masses of several times 10 6 M o and sizes 
of 100 pc (Blitz 1993). We adopt the nomenclature used 
in the review of Williams, Blitz & McKee (2000) to refer 
to the different observed molecular structures. GMCs 
are the sites of most of star formation activity in the 
Milky Way, and in particular of high mass stars, which 
are usually born in clusters within massive cores. In or- 
der to understand how massive cores form from the GMC 
complexes we must understand how fragmentation and 
condensation proceed within them. It is essential for this 
purpose to determine the physical properties of complete 
samples of massive cores within GMCs. 

Molecular line surveys, at millimeter and sub- 
millimeter wavelengths, have revealed the structure of 
GMCs to be highly inhomogeneous and clumpy (Blitz 
1993; Evans 1999). These surveys have shown that the 
mass spectra of clouds (Sanders et al. 1985; Solomon et 
al. 1987), clumps (Blitz 1993; Kramer et al. 1998; Lada 



1999; Williams, Blitz & McKee 2000) and total mass of 
embedded clusters (Lada & Lada 2003) are similar to 
one another. These mass spectra are notably different 
than the stellar mass spectrum: the Initial Mass Func- 
tion (IMF; Salpeter 1955. In particular, several works on 
molecular line mapping of GMCs, show that their mass 
spectra follow a power law with nearly the same expo- 
nent, x ~ 0.6, where dN/dlog M oc M~ x (e.g. Blitz 
1993; Williams, Blitz & McKee 2000). 

The development of large bolometer arrays during the 
last ten years has permitted to carry out extended mil- 
limeter and sub-millimeter continuum surveys, allowing 
the direct dust mass determination of clumps. Motte et 
al (1998) mapped the Ophiuchus cloud with the IRAM 
30-m telescope and found a distribution of clump masses 
similar to Salpeter's IMF. Several works (Johnstone et 
al. 2001; Beuther & Schilke 2004; Mookerjea et al. 2005; 
Reid & Wilson 2005; Johnstone et al. 2006) have re- 
ported mass spectra of dust cores with indexes similar 
to Salpeter's IMF, and different from those derived from 
molecular line studies. If the mass function of cores is 
similar to the Salpeter IMF, independent of the range of 
masses involved, then the star formation process within 
GMCs would be defined in the earliest stages as a result 
of cloud fragmentation. 
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In this paper we present a large scale 1.2 mm contin- 
uum study of the NGC 6334 GMC aimed to find and 
study a complete sample of massive cores. 



1.2. NGC 6334 and NGC 6357 

NGC 6334 is one of the nearest and most prominent 
sites of massive star formation, at a distance of only 1.7 
kpc (Neckel 1978). The central region of NGC 6334 con- 
sists of a ~ 10 pc long filament with seven sites of massive 
star formation. Within them there is a wide variety of 
activity associated with star formation, such as water 
masers, Hll regions (Rodriguez, Canto & Moran 1982; 
Carral et al. 2002), and molecular outflows. 

Early FIR (McBreen et al. 1979) and radio (Rodriguez, 
Canto & Moran 1982) surveys characterized the overall 
properties of the main sites of massive star formation, 
whereas NIR studies (Straw et al. 1989; Straw & Hy- 
land 1989) revealed the cluster forming nature of many 
of these sites. Kraemer & Jackson (1999) and Burton et 
al. (2000) review the different notations used by previous 
authors to identify bright sources. Sub-millimctric and 
millimetric results on NGC 6334 have been previously re- 
ported by Gezari (1982), Sandcll (2000) and McCutcheon 
et al. (2000), all of which focused in the northern por- 
tion of the main filament: sources I and I(N) which 
have been recently resolved into smaller cores (Hunter 
et al. 2006). In the case of Sandell (2000), the sources 
I and I(N) were redefined as the central peaks of each 
respective source. This after observations at 350 fim, 
450 /an, 380 fim, 1.1 mm and 1.3 mm showed inner 
structure for Gezari's cores. In Sandell (2000) sources 
I and I(N) have angular sizes of 10" x 8" and 11" x 8" 
respectively. These sizes are ~ 10 times smaller than our 
massive clumps ell and cl2 identified as I and I(N) re- 
spectively. For the present work, ell encloses sources I 
and I(NW) of Sandell's nomenclature, while cl2 encloses 
I(N), SM1,SM2,SM4 and SM5 as well as considerable 
extended emission in both cases. Hence, the comparison 
between that work and ours is not straighforward. 

Following the evolutionary sequence proposed by 
Beuther et al. (2006), NGC 6334 and NGC 6357 ap- 
pear to be in an intermediate phase of Massive Starless 
Clumps and Protoclusters. The latter region seems more 
evolved that the former since compact Hll regions are still 
prominent in NGC 6334 at radio wavelengths while the 
structure of NGC 6357 is more disrupted, suggesting that 
massive stars have already shaped the mother clouds. 
This is also supported by the presence of more infrared 
sources in NGC 6357. From this evidence, and the lack 
of significant amounts of cold material in between the 
two regions, we will consider NGC 6334 and NGC 6357 
as two independent regions, defining NGC 6334 as the 
southeastern portion of the map as shown in Figure 1. 

Even though NGC 6334 is one of the closest GMCs, 
it is still far compared to low-mass star-forming regions 
(e.g. Ophiucus, Orion B, Taurus). As a consequence, the 
detected clumps in NGC 6334 are larger and consider- 
ably more massive, and can be considered likely cluster- 
forming cores (Motte et al. 2003; Ward-Thompson et al. 
2006; Beuther et al. 2006) and we are unable to resolve 
their inner structure. 



2. OBSERVATIONS AND DATA REDUCTION 

The regions NGC 6334 and NGC 6357 were mapped 
using the 37 channel SEST Imaging Bolometer Array 
(SIMBA) in the fast mapping mode in three different 
epochs: July 2002, September 2002 and May 2003. The 
passband of the bolometer has an equivalent width of 
90 GHz and is centered at 250 GHz (1.2 mm). The 
half-power beamwidth of the instrument is 24" giv- 
ing a spatial resolution of 0.2 pc. Ninety five ob- 
serving blocks were taken towards the NGC 6334 and 
NGC 6357 regions, with typical extension of <~ 10' x 20', 
to sample a total area of ~ 2 square degrees between 
17H6 m 00 s - 36°40'00" and 17^28 m 00 s - 33°40'00". 
Skydip observations were done approximately every two 
hours to determine the zenith opacity at 250 GHz. Typ- 
ical opacities were r ss 0.2 with values ranging from 0.17 
to 0.4 in a few cases. We also checked pointing on r] Cari- 
nae every two hours and found a typical rms deviation of 
3-5" in azimuth and elevation. Every night we observed 
Uranus for flux calibration. 

The SIMBA data were reduced using the MOPSI pro- 
gram written by Robert Zylka after conversion by the 
simbaread program written at ESO. The SIMBA raw 
data consists of a time series for each of the 37 bolome- 
ters (channels) in the array. The time scries includes the 
counts per channel and sky position. The reduction pro- 
cedure first removes the brightest data spikes. Next a low 
order baseline in time is fit to the full observation file for 
each channel, and a zero order baseline is fit in azimuth 
for each channel. The data is then deconvolved by the 
time response function of each channel, as measured by 
the SEST staff. Gain elevation and extinction corrections 
are applied next. An iterative sky noise reduction proce- 
dure is then applied, where the counts of each channel are 
correlated with those of the other 36 channels, yielding 
a so-called flat field correction to calibrate the relative 
sensitivity of each pixel. A source image is finally pro- 
duced by averaging the flux of all channels as they pass 
through the same position on the sky. 

The sky noise reduction algorithm includes the flux 
coming from both the source and the sky simultaneously. 
If the source extends over several channels, this intro- 
duces spurious correlated flux which hampers the sky 
noise reduction procedure. To avoid this, a smoothed 
model of the source flux distribution on the sky is sub- 
tracted from the raw time series data. Thus, the sky 
noise reduction procedure can be repeated, finding a bet- 
ter source model with each iteration. 

The calibration was derived from maps of Uranus. The 
resulting multiplicative factor varied between 0.06-0.09 
Jy count -1 beam -1 . Finally, reduced images were com- 
bined with the MOPSI software to produce the final map 
(Figure 1). In order to reduce noise further, a Gaussian 
smoothing of 30" was applied to the final image and the 
map edges were removed. The typical rms noise of the 
final map is 25 mJy beam -1 . SIMBA observations usu- 
ally have an absolute flux uncertainty of ^20% (Faundcz 
et al. 2004). 

3. RESULTS 

3.1. Clump -finding Algorithms 

Different clump-finding algorithms have been used to 
study the substructure in molecular clouds. Among the 
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Fig. 1. — Grey scale image of the 1.2 mm emission towards NGC 6334 and NGC 6357. The polygonal line was used to exclude the noisy 
borders of the mapped region. The large square indicates the area defined as NGC 6334 in the present analysis. It encloses 182 of the 
347 clumps found using clfind2d. Also labeled are the three sub-regions chosen for statistical analysis: the central region NGC 6334a (see 
Figure 2); a larger extension NGC 6334b; and finally NGC 6334c which includes only the clumps outside of NGC 6334b. 



most used are Clumpfind (Williams, de Geus & Blitz 
1994) and Gaussclumps (Stutzki & Gusten 1990). Both 
algorithms have different biases, but find similar clump 
distributions (e.g. Schneider & Brooks 2004). We use the 
clumpfind algorithm because it makes no assumptions 
about inherent clump shapes. Clumpfind first finds the 
brightest emission peak in the image, then it descends to 
a lower contour level and finds all the image pixels above 
this level, associating them to the first peak (clump) if 
contiguous, or else defines one or more additional clumps. 
The spatial separation of clumps is defined along saddle 
points. We used a conservative lower detection threshold 
of 75 mJy (3cr) per beam. We found 347 clumps in the 
whole image, 182 of which are in the NGC 6334 region. 
Only two clumps are likely to be fictitious based on their 
small effective radii and location close to the borders, 
and only one of them is in the NGC6334 region. Thus, 
we will use the remaining 181 clumps in our analysis of 
the region NGC 6334 in this paper. Table 1 lists the 



clumps in our sample associated with previously known 
radio and IR sources. These sources are all located in 
the brightest part of the filament as seen in Figure 2. 

3.2. Clump Size Distribution 

Figure 3 shows a histogram of the clump size distri- 
bution. The effective radius, or size, is determined from 
the angular area encompassed by each clump (from the 
clfind2d output) assuming a distance to NGC 6334 of 
1.7 kpc. The sizes range from 0.1 to 1.0 pc, with a me- 
dian of 0.36 pc. These values are similar to those derived 
for clumps within GMCs (Blitz 1993; Williams, Blitz & 
McKee 2000; Pudritz 2002; Faiindez et al. 2004; Beuther 
et al. 2006). 

3.3. Mass Estimates 

Since the dust emission at 1.2 mm is most likely to be 
optically thin (e.g. Garay et al. 2002), the mass of each 
clump can be estimated from the observed flux density. 
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Fig. 2. — Peak position of the radio and FIR sources of Table 1 plotted over images of the SIMBA continuum at 250 GHz (left panel) 
and ATCA continuum at 1.6 GHz (right panel). 



TABLE 1 

Radio, FIR and Millimeter Sources in NGC 6334 



Radio 


Radio Position 


FIR 


FIR Position 


SIMBA 


SIMBA Position 


Ref 


Name 




Name 




Name 
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(J2000) 




(J2000) 




(J2000) 




G351. 20+0.70 


17 h 20 m 04.1 s -35°56T0" 






cl210 


17' l 19 m 58.0 s 


-35°55'56" 


(4),t 


unnamed a 




V 


17 h 19 m 57.4 3 -35°57'52" 


cl4 


17 h 19 m 56.7 s 


-35°57'56" 


(2),t 


A 


17 h 20 m 19.2 3 -35°54'45" 


IV 


17 ,l 20 m 20.7 3 -35°54'55" 


cl3 (cl7,cll5) 


17 h 20 m 19.8 s 


-35°54'44" 


(l),(2),t 


C 


17 h 20 m 32.6 s -35°51'24" 


III 


17 h 20 m 31.3 3 -35°51'48" 


clll 


17 h 20 m 34.2 s 


-35°51'32" 


(l),(2),t 


D 


17 h 20 m 44.3 s -35°49T8" 


11 


17 h 20 m 42.2 s -35°49'25" 


cll9 


17' l 20 m 42.8 s 


-35°49T6" 


(l),(2),t 


E 


17 h 20 m 50.9 s -35°46'06" 












(1) 


F 


17 h 20 m 53.4 s -35°47'02" 


I 


17 h 20 m 55.1 s -35°47'04" 


ell 


17' l 20 m 54.0 s 


-35°47' 00" 


(l),(2),t 






I(N) 


17 h 20 m 53.6 3 -35°45'27" 


c!2 


17 h 20 m 56.0 s 


-35°45T6" 


(3),t 



References. — (l)Rodriguez, Canto & Moran (1982),(2)McBrccn ct al. (1979), (3)Gezari (1982), (4)Moran ct al. (1990), f This work 
Note. — The FIR and Radio Equatorial coordinates have been processed from B1950 to J2000. 

a Thc radio counterpart for FIR source V is observed at 1.6 Ghz with similar peak intensity to G351.20+0.70, but was not identified earlier as 
an HII region but as part of a PDR. shell (Moran ct al. 1990; Jackson & Kracmcr 1999; Burton ct al. 2000) 

Td = 17 K 2 , we computed the mass of each clump found 



For an isothermal dust source, the total gas mass M g is 
related to the observed flux density at an optically 
thin frequency v as (c.f. Chini, Kriigel & Wargau 1987) 



S V D 2 



RdgHvB„(Td) 



(1) 



where k v is the dust mass absorption coefficient, D is the 
source distance, R c i g is the dust to gas mass ratio, and 
B v is the Planck function. In more convenient units the 
gas mass can be written as 



M g = 20.4 



S250GHZ 

Jy 

lcm 2 g~ 

«250GHz 



D 
kpc 



0.01 
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by clfind2d. We find that the clump masses in NGC 6334 
range from 3 to 6000 Mq, with a mean value of 170 A/ Q 
and a median value of 60 M Q . The total mass of the 
clumps in NGC 6334 is M tot - 5 x 10 4 M . 

3.4. Clump Mass Spectrum 

The Clump Mass Function (CMF), £ Q (M), is defined 
as the number of clumps per unit mass, 

dN _ AN 
dM ~ AM' 



(3) 



1 



Mr. 



(2) 



Using a dust mass absorption coefficient of K250 GHz = 
1 cm 2 g _1 (Ossenkopf & Hcnning 1994); a dust-to-gas ra- 
tio of 0.01; a distance of 1.7 kpc, and a dust temperature 



where AiV and AM are used to indicate observational es- 
timates. Some workers prefer to use a logarithmic CMF, 
£z , defined as the number of clumps per unit logarithmic 

2 We choose a value of 17 K because it lies in the typical range 
for cold dark clouds (Pudritz 2002). It is also also a factor of 
two smaller 34 K the average temperature of clouds with infrared 
counterparts (Faundez et al. 2004) making the comparison between 
both temperatures easier. 
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Fig. 3. — Clump Size Distribution (in pc) assuming a distance 
of 1.7 kpc to NGC 6334. The cross below the histogram indicates 
the median and standard deviation of the distribution. 
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(4) 



These two functions are related by the expression 
£ X (M) = (MlnlO)£ a (M). If £ a (M) has a power law de- 
pendence with mass, £ Q (M) oc M~ a , then the logarith- 
mic CMF should also have a power law dependence with 
mass, £ X (M) oc M~ x , where the power law exponents 
are related by x = a — 1. A Salpeter slope corresponds 
to x = 1.35 or a = 2.35. 

The CMF is usually estimated from a histogram of the 
derived clump masses assuming it has a power law form. 
The distinction between the observationally derived £ Q 
and £ x is made by Scalo (1998), Kroupa (2001) and Lar- 
son (2003) . Its implications in the type of binning used 
in the histogram is also mentioned in Klesscn & Burk- 
ert (2000). Figure 4 shows an histogram of the mass 
of the clumps within NGC 6334. The bin size A log M 
has a constant value of 0.5. The completeness limit is 
estimated to be ~ 30 M Q . A least squares linear fit 
to the AiV versus logM relationship gives a slope of 
-0.62 ±0.07. 

Variations of the CMF within NGC 6334.— Given the 
large number of clumps in our survey, we can assess 
possible changes in the CMF as a function of position 
within the cloud. We consider the three subrcgions 
within NGC 6334 defined in Figure 1 and construct a 
mass spectrum for each of them (see Figure 5) . The CMF 
in NGC 6334a, which encompasses the central filament 
containing the most massive clumps and where star for- 
mation is clearly taking place, is well fitted with a power 
law with an exponent of x = 0.11 ± 0.12, significantly 
shallower than the value determined for the whole sam- 
ple. The exponent steepens for NGC 6334b {x = 0.48), 
which covers a much larger region than the central fila- 
ment. For NGC 6334c, which excludes the main filament, 
the slope is even steeper, x = 0.82. Thus, the slope of 
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Fig. 4. — Mass histogram of the clumps within NGC 6334. AN 
is the number of clumps in bins of constant AlogM (=0.5). The 
shaded region delimits the 3<x detection limit and our estimated 
completeness limit. This region is excluded from the least squares 
fit. 

the CMF in NGC 6334 depends on the location of the 
clumps within the cloud. 

Out of the total mass of 5 x 10 4 Mq in 181 dense mas- 
sive clumps in NGC 6334, 3.4 x 10 4 Mq are contained 
within the NGC 6334a region and 1.6 x 1O 4 M are out- 
side. The slope of the CMF of the whole region is dom- 
inated by the relatively low mass clumps in the outer 
region since they dominate by number. The bulk of the 
cloud mass is located in a few inner clumps which do 
not affect significantly the exponent of the derived CMF. 
This could shed light on the effects of the star formation 
activity, mass segregation, or coalescence in the clump 
mass spectrum within GMCs. We know that it is more 
likely to miss low mass clumps due to confusion within 
the region NGC 6334a than in the outer regions under 
study. This bias naturally flattens the slope of the mass 
spectrum in the inner region. However, we find in §3.6.1 
that clumps are not only concentrated by mass towards 
the center, but also by number, hinting that this result is 
partly real. Higher angular resolution observations will 
be needed to settle this issue. 

3.5. Possible Uncertainties in the CMF 

In the above analysis we assumed a single temperature 
for the whole ensemble of clumps. This is clearly an ap- 
proximation; the clumps themselves are not isothermal 
and the temperature is likely to be different from clump 
to clump. In addition, we assumed that all the detected 
1.2-mm emission is due to dust thermal emission. It is 
possible, however, that some of the 1.2-mm emission is 
due to free-free emission from ionized gas. In what fol- 
lows we assess these two assumptions and quantify their 
effects on the derived CMF. 

3.5.1. Temperature 

Assuming that all clumps within a GMC are isothermal 
and have the same temperature is clearly a rough approx- 
imation. In particular, clumps with already formed stars 
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Fig. 5. — Mass histogram of clumps within selected regions of NGC 6334 (see Figure 1). Top panel: NGC 6334a. x = 0.11 Middle panel: 
NGC 6334b. x = 0.48 Bottom panel: NGC 6334c. x = 0.82 
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Fig. 6. — Three color image of the MSX emission towards NGC 6334 (blue: Band A, 8.28/im; green: Band C, 12.13/^m; red: Band E, 
21.3/xm). The countours represent the SIMBA 1.2-mm emission. Green contours indicate clumps that are associated with extended or 
point-like infrared counterparts. 
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Fig. 7. — ATCA radio continuum emission at 1.6 GHz overlaycd 
with contours of the 1.2-mm emission. Again, the red contours 
indicate cold clumps having no embedded MSX sources and the 
green contours indicate warm clumps with embedded MSX sources. 

are expected to be warmer than clumps with no signs of 
embedded objects. Beuther & Schilke (2004) have argued 
that this assumption introduces an uncertainty in the de- 
rived slope of the mass spectrum. If higher temperatures 
are adopted for the more massive clumps their derived 
masses will decrease, whereas if lower temperatures are 
adopted for the less massive clumps their derived masses 
will increase: the change would steepen the slope of the 
spectrum. 

To study the possible effects of temperature differ- 
ences, we use MSX mid-infrared and/or ATCA cm- 
continuum observations to determine the presence of em- 
bedded heating sources which arc likely to be respon- 
sible for temperature differences between clumps. Fig- 
ures 6 and 7 show, respectively, images of the MSX and 
ATCA emission overlayed with contours of the SIMBA 
1.2-mm continuum emission. We identified by visual 
inspection clumps associated with extended or point- 
like infrared counterparts (Figure 6), or clumps associ- 
ated with significant radio continuum (free-free) emission 
(Figure 7). These clumps are likely to have embedded 
sources and are indicated with green contours in both 
images. Clumps in red contours appear to be free of 
embedded infrared or radio sources. 

For the purposes of constructing the CMF, the power- 
law index does not depend on the temperature chosen for 
the whole ensemble of clumps, and the choice of higher 
temperatures only displaces the histogram to the left. 
Any slight variation of x is due to the intrinsic prob- 
lem of binning the data. Nevertheless, the choice of two 
different temperature can change the shape of the his- 
togram, but essentially depending on the ratio between 
the temperatures chosen rather than the values them- 
selves. 

We identified 16 out of 181 cores that appear to be 
warmer than the rest. For these clumps we adopt a tem- 
perature of 34 K, the average temperature of massive 
clumps with embedded IRAS sources as determined by 
Faundez et al. (2004). In particular, the temperatures 



assigned to the cores associated with NGC 6334 I and 
NGC 6334 I(N), of 34 K and 17 K respectively, are in 
accordance with the values given by Gezari (1982). How- 
ever, Sandell (2000) argues that source I is much hotter 
(T d ss 100 K) while his estimate for I(N) is 30 K, just 
a factor of two larger than our estimate. His estimate 
of the mass of NGC 6334 I of 200 M & in contrast to 
the 1800 Mq in the present work is not solely explained 
by the temperature difference but also by the integrated 
flux. With a resolution of 6" at 800 /im, his estimate for 
the size of source I is 10" x 8", implying an effective radius 
~ 8 times smaller than our estimate. Similarly, the defi- 
nition of Sandell (2000) of NGC 6334 1(N) corresponds to 
the central peak of the larger clump detected by Gezari 
(1982). Sandell (2000) obtains a mass of 400 M for 
source I(N) after assuming a temperature of 30 K (see 
Sandell 2000 and references therein) . Adding up the con- 
tributions from the different cores resolved within I(N) 
and including the surrounding cloud, Sandell (2000) finds 
that the mass of the whole I(N) region is ~ 2700 M Q This 
value is in agreement with other results found in the lit- 
erature and is consistent with our result after taking into 
account that our choice for the temperature is 17 K. We 
remark that estimating an exact value of the tempera- 
ture for each of the clumps of the cloud is not relevant for 
the statistical analysis we carry out in the present work. 

The newly adopted temperatures imply a repositioning 
of ~10% of the clumps in the mass histogram. The new 
mass histogram, made assuming a two temperature cloud 
ensemble, is shown in Figure 8. A linear regression yields 
a best fit value of x — 0.88±0.13 (dotted line), but a \ 2 fit 
with Poissonian error bars yields x = 0.62 ± 0.08 (solid 
line). Thus, even though the warmer clumps are also 
preferentially the most massive, they do not concentrate 
solely on the most massive bin and do not affect the 
derived slope of the mass spectrum significantly. 

The repositioning of clumps in the histogram could 
have a more dramatic effect when the number of 
clumps is considerably smaller. For example, the area 
NGC 6334a includes 14 of the 16 warmer clumps in 
NGC 6334 and a total of only 40 clumps (figure 5). 
When using two temperatures in NGC 6334a, the shape 
of the histogram indeed steepens, but the large relative 
errors yields a slope x = 0.35 ± 0.21. Excluding the 
last bin, containing only one object, we obtain a slope 
of x = 0.09 ± 0.15. In both cases the slope remains con- 
sistent with the value x ~ 0.1 from a single temperature 
clump mass distribution. 

3.5.2. Radio Emission Contribution 

Free-free emission is the main continuum contributor 
at radio frequencies. At frequencies of hundreds of GHz, 
the contribution from bremsstrahlung emission is usually 
neglected and all the detected emission is assumed to be 
due to dust. Since massive stars are being formed within 
NGC 6334, it is plausible that the 1.2-mm emission is 
not completely due to cool dust emission and that the 
ionized gas excited by these stars contribute an impor- 
tant amount. Brooks et al. (2005) studied the Keyhole 
nebula and found that the 1.2 mm emission towards the 
Hll region Car-II is strongly correlated with the 4.8 GHz 
continuum emission and that there is a lack of molecular- 
line emission. They concluded that the 1.2-mm flux from 
the components of Car-II arise from free-free emission as- 
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Fig. 8. — Histogram of mass distribution in NGC 6334 for a 
two-temperature clump ensemble: 17 K and 34 K. 



TABLE 2 

Gaunt Factors for Free-Free Emission 



Frequency 


T=7000 K 


8000 K 


9000 K 


10000 K 


1.6 GHz 


5.40327 


5.5134 


5.61051 


5.69735 


250 GHz 


2.79429 


2.8871 


2.9702 


3.04545 



sociated with ionized gas and not from cool dust emission 
associated with molecular gas. 

Here we make a correction to the observed emission 
at 250 GHz by free-free contamination by estimating the 
expected ionized gas flux density at 250 GHz from the 
observed flux density at 1.6 GHz. Assuming that the 
free-free emission is optically thin at both frequencies, 
the ratio of the emissivities is proportional to the ratio 
of the e~ hv l hT,1 gff(y,T e ) factors (Rybicki & Lightman 
1979), where gjf is the Gaunt factor. The exponential 
is essentially 1 at both frequencies. We did not use the 
usual radio approximation for the Gaunt factors given 
by Altenhoff ct al. (1961), but computed them more pre- 
cisely using quantum mechanical calculations following 
the work of Menzel & Pekeris (Menzel & Pekeris 1935; 
Sommerfcld 1953). Table 2 lists the calculated gaunt fac- 
tors averaged over a Maxwell-Boltzmann distribution of 
velocities with temperatures between 7 x 10 3 and 10 4 K, 
typical for Hll regions in massive star forming regions 
(Beckert et al. 2000). At 1.6 GHz the computed Gaunt 
factors are only w 0.3% lower than the usual Altenhoff 
ct al. approximation. On the other hand, the calculated 
Gaunt factors at 250 GHz are typically 15% smaller than 
the values obtained from the Altenhoff et al. approxima- 
tion. 

Table 3 summarizes the values of the estimated free- 
free emission at 250 GHz from selected clumps associated 
with HII regions (see Figure 2). We subtracted this 
flux from the measured 1.2-mm flux density to estimate 
the actual contribution from dust, rederiving the mass of 
each clump. The high mass bins are the most affected, 



but the best fit exponent [x 
with the previous value. 



3.5.3. 



0.88 ± 0.13) is consistent 



Binning 

Given a distribution of clump masses, the binning pro- 
cess invariably loses information (Rosolowsky 2005). The 
mass histogram can be interpreted as the derivative of a 
cumulative number function Af(M) which counts clumps 
with mass greater than M, 



Jo £,a(m)dm 



(5) 



Some authors argue in favor of using the cumulative num- 
ber function Af(M) (Johnstone et al. 2000b, 2001; Ker- 
ton et al. 2001; Tothill et al. 2002) to avoid the loss of 
information. If the cumulative number function has a 
power-law dependance with mass, M(M) cx M~ 7 , then 
the CMF is also a power law £ Q (M) cx M~ 7_1 . From 
the definitions in section 3.4, it is clear that 7 — x 

However, if there is an upper mass limit in the distri- 
bution of cores, then the cumulative mass function is not 
a power-law , showing considerable curvature at the high 
mass end. Af(M) can be approximated by a power-law 
with index x = a — 1 only at masses M <C M max (or 
when M max —> 00). For finite upper mass limit, 



Af(M) = CiM- x + C* 2 



(6) 



where C'2 becomes unimportant for small masses (see 
Rosolowsky 2005; Reid & Wilson 2006b; Li et al. 2006). 
A power-law form, A/"(M) cx M~ x 7 approximates the 
true CMF asymptotically towards low masses. Thus, 
a power law fit to the CMF must be applied within a 
range of masses that avoids both the incomplete low- 
mass range and the cutoff high-mass range. When fitting 
a power-law function to an observed cumulative mass 
function in a mass range M > M max /2 the slope in- 
variably increases, explaining the apparent Salpcter-like 
slopes found in previous studies. Furthermore, the slope 
thus obtained is strongly dependent on the break-point 
chosen to fit the high mass end of the cumulative mass 
function and does not reflect the underlying differential 
clump mass distribution (£ Q or £ x ). 

For sample sizes of ~ 70 clumps or fewer, binning 
becomes an important factor in fitting a power law to 
differential mass functions and the use of a cumulative 
mass function is preferred. Johnstone et al. (2000b, 2001, 
2006) and Reid & Wilson (2005, 2006a) use cumulative 
mass distribution functions to avoid this problem. Our 
sample (181 clumps) is large enough to analyse the data 
using either the cumulative or the differential mass func- 
tions. Figure 9 plots the normalized cumulative number 
function of clumps in NGC6334. A single power-law fit, 
which we have shown does not represent the true under- 
lying clump mass function, gives Af(M) cx M~ 0,87 . The 
dashed line in Figure 9 shows that the slope x = 0.62 
derived from the histogram (Figure 4) represents a good 
asymptote to the CMF, as predicted by the theory (Eq. 
6). We also show in Figure 9 a fit to the top 10% of 
the cloud mass (log(M/M Q ) > 2.8, indicated by an ar- 
row) which yields a slope x = 1.19 (dotted line) much 
closer to Salpeter's value. The best fit exponent changes 
to x = 1.37 for logM > 2.9, and to x = 1.47 for 
logM > 3.0. Even though these slopes are consistent 
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TABLE 3 

Radio an MM fluxes [Jy] for Selected Clumps in Figure 7 



Frequency 


ell 


cll9 


cll69 


clll 


cl27 


cl7 


cl210 


cl4 


ATCA at 1.6 GHz (Free-Free) 
Expected Free-Free at 250 GHz 
SIMBA at 250 GHz 


0.59 
0.31 
68.38 


22.77 
12.16 
26.78 


0.45 
0.24 
1.13 


9.42 
5.03 
27.12 


5.28 
2.82 
13.69 


12.61 
6.73 
26.72 


3.49 
1.86 
1.34 


2.89 
1.54 
59.31 



1 .00 



0.10 I 



0.01 ■ 




1 .0 



2.0 2.5 

log(M/M ) 



3.0 



3.5 



Fig. 9. — Normalized cumulative mass function Af(M). The grey area corresponds to the region affected by completeness uncertainties 
as shown in Figure 4. The arrow delimits the top 10% of the CMF where a power-law fit to the CMF (dotted line) yields an apparent 
Salpcter-like slope. 



with Salpeter's IMF, they are an artifact of having an 
upper clump mass limit in the sample and do not reflect 
the true clump mass function 

In summary, one must be aware of introducing biases 
when using the high mass range to fit power laws or when 
fitting broken power laws to the cumulative mass func- 
tion. These problems are minimized (but still present) 
with samples larger than 100 clumps. 

3.6. Spatial Distribution of Clumps 

In order to understand the process of fragmentation, 
we need to explain how masses are distributed in clumps 
and how they are positioned in space. A complete theory 
of star formation must not only reproduce the mass func- 
tion, but it must explain it in all its physical implications 
including how clumps, cores and stars are distributed 
spatially during the evolution of the GMC (Bonncl ct 
al. 2006). The exhaustive study of the spatial distribu- 
tion of young stars in the Taurus region by Gomez et 
al. (1993) was extended by other authors and compiled 
by Larson (1995). Here we undertake a study the de- 
gree of clustering of clumps in NGC 6334 using a similar 
approach. 

3.6.1. Clustering and Segregation of Clumps 



The number density of clumps can give us insight 
about the actual state of fragmentation and how the 
clumps are distributed spatially independent of their 
mass. We study the number density of clumps by means 
of the Simple Grid 3 and the Kernel Methods 4 . Both 
methods require a free parameter which determines the 
"resolution" of the number density estimator: the bin- 
ning length D in the grid technique, and the smoothing 
length h in the Kernel Method, where a kernel K is de- 
fined at each pixel of the map in Figure 1 by 



2tt 



with the kernel density estimator D defined by 



1 ™ 

D{x, y) = ^2 K ( x ' x ^ y> Vi) 



(8) 



where x and y are measured in pc ignoring the sky cur- 
vature. The grid bin / was taken to be 2 pc while h was 

3 The Grid Method consists in binning the two-dimensional 
space with squares of side I and then dividing the numbers of 
sources lying within each square by the area of it I 2 , to obtain 
number density in units of length -2 

4 The Kernel Method (Silverman 1986) uses a kernel function K 
offering the advantage that the density distribution is smoothed. In 
each point (x,y) of (a,<5), the kernel density estimator determines 
the density due to the contributions of all n data points. 
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Fig. 10. — Normalized Nearest Neighbor Distribution after di- 
viding by 181, the number of neighbor pairs (solid line) and the ex- 
pected distribution for random (Poisson) distribution of the same 
number of objects as the original sample over an identical area 
(dot-dashed line). The medians of the distributions are indicated 
by dotted lines. 



chosen in order to smooth the distribution over an area 
irh 2 similar in size to the area over which the simple grid 
technique smoothed the data points (I 2 ). 

We applied this analysis to the subregion NGC 6334b 
which covers and area of ~ 330 pc~ 2 . We find that the 
probability that the clumps are distributed at random 
within this subregion is ~ 10 , which argues in favor of 
clustered fragmentation at scales between 0.1 and 10 pc. 
Thus, we conclude that there is spatial segregation in 
clump number in addition to clump mass. 

3.6.2. The Nearest- Neighbor Distribution in NGC 6334 

We calculate the nearest-neighbor distribution (i.e. the 
frequency distribution of the linear distance to the near- 
est neighbor of each clump) for the clumps in NGC 6334. 
We neglect those clumps located too close to the edge of 
the mapping area but do not exclude them from the total 
sample since they can be the nearest neighbor for an in- 
ner clump. We binned the nearest neighbor distances in 
intervals of 0.2 pc to construct the histograms shown in 
Figure 10. The nearest neighbor distribution is strongly 
skewed to small separations and very different from the 
distribution expected from random positions at the same 
mean density. The differential probability of observing 
at least one event in the interval [r\ , r 2 ] which defines a 
ring surrounding a central source is (Gomez ct al. 1993) 



AP{n,r 2 ,v) = e~^ - e - nr > r * 



(9) 



where r\ = N/A is the average density, in this case given 
by n = 0.11 clumps pc~ 2 . Figure 10 compares the Pois- 
son distribution histogram (dashed line) with the nor- 
malized nearest-neighbor distribution for clumps. The 
median for the random distribution (1.7 pc) is a fac- 
tor of <~ 2.5 greater than that derived from the actual 
clump distribution (dotted vertical lines in Figure 10). 
We would need to increase the average density rj by a 
factor of 5 to produce a random distribution with the 
same median separation as the actual clump distribu- 
tion. Thus, we conclude that the median separation of 
the clump distribution is significantly smaller than the 
expected median separation of a random distribution. 




Fig. 11. — Average Surface Density as a function of separation 
in degrees (Equation 10). The vertical dotted lines encompass the 
limits of validity for the fit: the beamsize logfi = —2.7 (~ 0.2 pc), 
and the sampling bias (the width of subregion NGC 6334b) at 
log 6» = -0.5 or 10 pc. 

3.7. The Average Angular Surface Density of Clumps in 

NGC 6334 

Most of the newly formed stars in nearby regions of 
star formation are located in groups or clusters. The 
degree of clustering of prc-main sequence stars can be 
obtained by measuring their surface density as a function 
of angular distance, £(0), from each star. This surface 
density can be fit by power-laws (Gomez et al. 1993); 
Larson (1995) found that there is a characteristic spatial 
scale (~ 0.04 pc) where the surface density changes slope. 
This scale could mark a transition between the regime 
of cores within molecular clouds and protostars within 
cores. 

We applied the approach followed by Larson (1995) 
to the clumps in NGC 6334, which have typical sizes 
~ 0.1 pc and separations of <~ 1 pc. We compute £(#) 
by taking each clump k and dividing the surrounding 
area of the sky into a set of annuli of radius 6^ (with 
6i = 1.5#j_i), and counting the number of companion 
clumps Nk{0i) in each annulus (Kitsionas et al. 1998). 
Then. 



£(00- ' 



{6j + Oj 
2 







(10) 



The results are plotted in Figure 11. The surface den- 
sity is poorly fitted by a single power-law, but the fit is 
much improved when using a broken-power law. How- 
ever, this reflects a sampling problem at large separa- 
tions, since the area observed around NGC 6334 is not 
square, Its narrowest part has a width of ~ 10 pc, just 
where the power-law apparently breaks. The reliable 
portion of the power-law has a slope of —0.62 



£ c oc 6 



-0.62 



0.6 pc < 6D < 10 pc 



(11) 



which is remarkably similar to the one found by Larson 
(1995) for separations larger than 0.04 pc calculated from 
young stars in the Taurus-Auriga region: 

-0.62 



XL cx 



0.04 pc < 6D < 2.5 pc. (12) 



These two results cover two different separation ranges 
in regions with widely different physical properties, yet 
the surface density distribution is the same. If cluster- 
ing above the characteristic length of 0.04 pc maintains 
a self-similar behavior up to 10 pc, then the grouping 
of stars at large separations (between subclusters for ex- 
ample) could be determined from early stages of cloud 
fragmentation before stars are formed. 
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4. DISCUSSION 

The CMF and its relation to the IMF. — The clump mass 
function is, on one hand tied to the formation and evo- 
lution of their parent molecular clouds, and on the other 
to the formation of their daughter stars. Provided that 
the mechanism of clump fragmentation and collapse to 
form stars is universal, then the IMF should be a di- 
rect consequence of the CMF. However, since a single 
clump generally forms multiple stars, we can't expect 
the distribution of clump masses and stellar masses to 
have identical functional form. 

We find that the mass spectrum of clumps in NGC 6334 
has a power-law dependence with mass with an index 
x = 0.62 (dN/dlogM oc M~ - 62 ). This value is similar 
to values derived from dust continuum observations for 
other massive star forming regions (Kerton et al. 2001; 
Tothill et al. 2002; Mookerjea et al. 2005). It is also sim- 
ilar to values derived from molecular line observations, 
usually isotopomeric lines of CO, for clouds with simi- 
lar total mass. For instance, Nozawa et al. (1991) finds 
a = 1.7 in Ophiucus North (mass range of 4-250 M©), 
Stutzki & Gusten (1990) finds the same exponent in 
M17SW (mass range of 10-3200 M Q ), and Kramer et 
al. (1998) reports a = 1.8 for NGC 7538 (mass range of 
50-3.9 x 10 3 Mq). These results indicate that the clumps 
in all these massive star forming regions are not the di- 
rect progenitors of individual stars. The fact that the 
massive clumps follow a power-law mass function over a 
wide range in mass is remarkable. The similarity of the 
"high mass" CMF slope to the Mass Function of GMC's 
in our Galaxy (Sanders et al. 1985; Solomon et al. 1987; 
Pudritz 2002) seems to support a hierarchical or fractal 
model of the distribution of gas in the Milky Way, where 
fragmentation and mass distribution can be interpreted 
as scale free. 

In order to build a complete theory of star formation we 
must understand the process of clustered star formation 
in clumps satisfying the observed CMF (dN/dlogM oc 
M~ - 6 ) and leading to the observed IMF (dN/d\og M cx 
M -1 - 35 ). This has stimulated studies about the process 
of fragmentation of GMCs and clumps both theoretically 
and observationally. Hydrodynamic and magnetohydro- 
dynamic numerical calculations suggest that turbulence 
might play a major role in the clump fragmentation pro- 
cess. For example Bonncll & Davies (1998) suggest that 
stars are formed through the competitive accretion of gas 
onto proto-cores within molecular clouds. 

Recent sub-millimeter and millimeter continuum ob- 
servations of low-mass star forming regions (e.g. Ophi- 
ucus, Serpens, Orion B) (Motte et al. 1998; Testi & 
Sargent 1998; Johnstone et al. 2000b, 2001; Beuthcr & 
Schilke 2004) have revealed a CMF with a slope similar to 
that to the IMF. The different values obtained from the 
molecular line observations described above and these 
continuum observations probably reflects the different 
mass ranges sampled. In fact, dust continuum emission 
observations toward more distant massive star forming 
regions have revealed CMF slopes consistent with a value 
of -0.6 (Kerton et al. 2001; Tothill et al. 2002; Mooker- 
jea et al. 2005; this work), similar to that obtained from 
molecular line observations. From observations of a sam- 
ple of massive clumps toward the massive star forming 
region NGC 7538, Reid & Wilson (2005) concluded that 



a Salpeter-like mass function is already established at 
the earliest stages of star formation. This is a surpris- 
ing result since many of their clumps, with masses of 
10 2 — 1O 3 M , are still likely to be undergoing the pro- 
cess of fragmentation and can not be direct progenitors 
of individual stars. Even if the structure in ISM were 
fractal, self-similarity must break on small scales, where 
star formation is taking place. We note, however, that 
the result of Reid & Wilson (2005) comes from a fit to 
the high mass end of the cumulative mass function. We 
showed in section 3.5.3 that these fits are biased towards 
larger exponents. In addition to molecular line and dust 
continuum surveys, extinction maps can be used to map 
the dense molecular cores in star forming regions. Us- 
ing this technique, Alves et al. (2007) have found a mass 
spectrum for cores in the Pipe Nubula that is surpris- 
ingly similar to the IMF. Their CMF displaced to higher 
masses with respect to the stellar IMF only by a factor 
of 3, suggesting a one-to-one mapping from cores to stars 
with a star formation efficiency of <~ 30%. In their case, 
the cores detected span a range of masses of 1-10 M Q ap- 
proximately This clearly indicates a much smaller scale 
than the clumps in our work whether determined by the 
tracer or just the spatial resolution of the observations. It 
is not clear if the scale of the observations determines the 
observed fragmentation conditions i.e. the distribution 
of masses or if the IMF is a result of the complex evo- 
lution of the accreting cores and their interplay with the 
harboring molecular cloud and other companions. The 
low star forming activity in the Pipe Nebula has led ? 
to suggest that this CMF can be determined at early 
evolutionary stages. 

Assuming that at scales of 0.4 pc and masses from 
1OOM to 5OOOM the Blitz slope (= 0.6) is valid, the 
question that arises is: What must happen to change 
the slope from Blitz-like to Salpeter-like? In the context 
of gravitational opacity limited fragmentation, gravita- 
tional collapse starts from density inhomogeneities and 
proceeds with cooling, which in turn produces smaller 
Jeans masses in the colder regions, favoring gravita- 
tional fragmentation on small scales. In a strictly self- 
similar regime, fragmentation should occur maintaining 
the same slope in the mass spectrum at any scale. But 
at some point the gas cores will not be supported by 
thermal or non-thermal motions, collapse will occur and 
stars will form, halting the clump fragmentation, while 
larger and less dense clumps will continue fragmenting. 
Eventually large clump masses will be depleted, and the 
number of small mass clumps will increase steepening 
the slope to eventually reach a Salpeter value. The exact 
form in which this happens is probably a combination of 
many of the mechanisms proposed. 

Will these clumps form stars? — The discussion 
sumes that all clumps fragment to form stars in order to 
obtain the IMF from a CMF. However, this might not 
be the case, in particular for some of the least massive 
clumps further away from the cloud center. Are these 
clumps gravitationally bound? Are they hence likely to 
collapse or are they only transient structures or overden- 
sities triggered by turbulent compressive shocks? These 
questions have yet to be answered,either observationally 
or theoretically. 
Numerical simulations commonly report that many 
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of the lower-mass cores formed are not gravitationally 
bound (Klessen 2001; Clark & Bonnell 2005; Tilley & 
Pudritz 2005) . Furthermore, the mass spectra can be un- 
derstood as due to purely hydrodynamical effects with- 
out gravity (e.g. Clark & Bonnell 2006). As Padoan 
& Nordlund (2002) remark, the mass spectra resulting 
from turbulent fragmentation is different from the one 
that considers collapsing or unstable cores. Only the lat- 
ter form stars. These calculations show that many of the 
clumps could be transient structures; indeed a significant 
fraction of the cores end up re-expanding rather than col- 
lapsing (Vazquez-Semadeni et al. 2005; Nakamura & Li 
2005). This implies that fragmentation is not sufficient 
to trigger star formation. If supersonic turbulence gen- 
erates the initial density enhancements from which cores 
develop, then these cores might not necessarily approach 
hydrostatic equilibrium at any point in their evolution 
(Ballesteros-Paredes et al. 1999). 

Can we distinguish observationally between bound and 
transient clumps? To estimate if clumps are likely or 
not to collapse, we need to know more than its mass. 
A clump's column density is important as a diagnostic 
of whether the clump is likely to collapse (Pudritz 2002) . 
Values higher than TV ~ 10 22 cm~ 2 are observed for cores 
with embedded sources. The virial parameter, M v ; r /M, 
where M is a mass derived from column density and M v ; r 
is the virial mass derived from the cloud radius and veloc- 
ity dispersion(Bertoldi & McKee 1992), is known to have 
values close to 1.0 in star forming clouds (Onishi et al. 
1996; Yonekura et al. 2005). However, given the observa- 
tional uncertainties it is not yet clear if this parameter is 
a good diagnostic for star forming vs. transient clumps. 

Preferred spatial scale. — As Larson (1995) found in Tau- 
rus, there might exist a preferred scale of star formation 
at which clustering changes. This length can be related 
to the Jeans length and the Jeans mass. The surface 
number density of the clumps as a function of separation 
can reveal a characteristic spatial scale, marking a tran- 
sition between clumps and cores. At separations larger 
than 0.6 pc and smaller than 10 pc, we find the same 
slope for the power-law fit to the surface density of com- 
panions as the protostars in Taurus do at large separa- 
tions: E cx 9~ - 62 . This suggests that large separations 
in stellar systems are determined by the position of their 
progenitor clumps. 

5. CONCLUSIONS 

We made observations of the 1.2-mm dust emission 
toward the Giant Molecular Cloud NGC 6334 using the 
SIMBA bolometer at the SEST. The main results and 
conclusions presented in this paper are summarized as 
follows. 

We find 181 clumps, distributed in an area of <~ 2.0 
square degrees centered on the main filament, which har- 
bors most of the star-forming activity. The clumps range 
in size from 0.1 to 0.9 pc, with a median of 0.35 pc. This 
range is similar to that found by Faiindez et al. (2004) 
and Plume et al. (1997) for clumps in different high mass 



star forming regions. The clump masses, assuming they 
are isothermal, range from 3 M Q to 6 x 10 3 M Q , with a 
completeness limit of ~ 30 M Q (assuming T<j = 17 K). 

The Clump Mass Function (CMF) is well fit with a 
power-law dependence with mass with an index x = 
0.62 ± 0.07 (dN/dlogM cx M~ x ) in the mass range be- 
tween 30 M to 6 x 10 3 M Q . The slope differs from the 
stellar IMF slope, indicating that clumps are not direct 
progenitors of stars. Therefore other processes besides 
fragmentation must be important in setting up the IMF 
from the CMF. 

We assessed possible effects on the derived slope of the 
CMF due to changes on the temperature assumptions 
and due to the contribution of free-free emission from 
ionized gas to the 1.2-mm emission. Although < 10% 
of the clumps arc likely to be significantly warmer than 
17 K and are associated with regions of ionized gas, the 
correction for temperature and free-free emission has lit- 
tle effects on the derived slope. 

We investigated possible differences in the value of x 
among different sub-regions of NGC 6334. We find that 
the slope is significantly shallower toward the central fila- 
ment [x w 0.1), which contains the most massive clumps 
and represents the minimum of the gravitational poten- 
tial in the GMC. As we cover more extended regions, 
with clumps not actively forming stars, the slope steep- 
ens (« 0.5, 0.8), revealing that the bulk of the clumps 
arc located in the outer areas of the molecular cloud and 
that these low-mass clumps predominantly determine the 
shape of the mass function. 

We caution about the power-law fitting procedures to 
the mass function. The differential CMF is sensitive to 
bin size and to low-number statistics in the last bin (high- 
mass end) as well as to completeness limits in the low- 
mass end. On the other hand, fitting a power-law to the 
high mass end of the cumulative CMF is incorrect due to 
its high mass cutoff. Both the low-mass and high-mass 
ends of the cumulative CMF must be avoided in fitting 
power-laws. 

The spatial analysis performed on the two-dimensional 
distribution of clumps reveals that they are not dis- 
tributed randomly. They are concentrated toward the 
center of the filament, indicating not only a segregation 
in mass but also a segregation in number which could 
suggest a possible coalescence of massive clumps towards 
the gravitational potential minimum. In addition, we 
study the surface density of companions as a function of 
separation. This is well fit by a power-law with a similar 
exponent to the one found for proto-stars in Taurus at 
large angular separations. This suggests that the posi- 
tion of stars in clusters is determined in the fragmenta- 
tion and star formation stage rather than after dynamical 
relaxation. 

We thank Simon Casassus and Dieter Nurnberger for 
helpful comments. D.J.M., D.M., G.G., and D.R. grate- 
fully acknowledge support from the Chilean Centro dc 
AstrofMca FONDAP No. 15010003. 
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TABLE 4 
Clump Properties in NGC 6334 
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17' , 17 m 58.5 s 


-36° 9' 6.1" 


0.24 


2.00 


0.42 


1.23X10 2 


1.23X10 2 


7.21x10 s 


NGC 


63 34(6334c) 
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TABLE 4 
Clump Properties in NGC 6334 



cl 132 


17 fe 18 m 27.2 s 


-36°25'39.7" 


0.24 


1.28 


0.36 


7.90X10 1 


7.90X10 1 


7.34x10 s 


NGC 


63 34(6334c) 


cl 133 


I7 h 20 m 7.9 s 


-35°58'45.8" 


0.23 


0.92 


0.30 


5.66X10 1 


5.66X10 1 


9.09x10 s 


NGC 


6334(6334a)(6334b) 


cl 134 


17 h 22 m 3.7 s 


-35°24'42.1" 


0.23 


1.99 


0.45 


1.23X10 2 


1.23X10 2 


5.85x10 s 


NGC 


6334(6334b) 


cl 135 


17 h 22 "i i.8 s 


-35°26T8.2" 


0.23 


1.09 


0.33 


6.71X10 1 


6.71X10 1 


8.09x10 s 


NGC 


633A(6334b) 


cl 136 


17 h 21 m 32.4 s 


-35°31'31.8" 


0.23 


1.49 


0.42 


9.19X10 1 


9.19X10 1 


5.37x10 s 


NGC 


6334(6334b) 


cl 138 


17 h 19 m 28.9 s 


-36° 9'57.6" 


0.23 


1.56 


0.42 


9.64X10 1 


9.64X10 1 


5.64x10 s 


NGC 


6334(6334b) 


cl 140 


17 fc 21 m 33.1 s 


-35°32' 3.8" 


0.22 


0.74 


0.27 


4.55X10 1 


4.55X10 1 


l.OOxlO 4 


NGC 


6334(6334b) 


cl 141 


17 h 18 m 26.8 s 


-36°14T1.8" 


0.22 


1.24 


0.36 


7.68X10 1 


7.68X10 1 


7.14x10 s 


NGC 


6334(6334c) 


cl 142 


17 h 21 m 6.4 s 


-35°47'33.0" 


0.22 


0.56 


0.24 


3.47X10 1 


3.47X10 1 


1.09X10 4 


NGC 


6334(6334a) (6334b) 


cl 143 


17 h 17 m 58.4 s 


-36° 10' 2.3" 


0.22 


1.30 


0.36 


8.06X10 1 


8.06X10 1 


7.49x10 s 


NGC 


63 34(6334c) 


cl 145 


17 fc 21 m 11.7 s 


-35°50' 4.6" 


0.22 


0.86 


0.30 


5.34X10 1 


5.34X10 1 


8.57x10 s 


NGC 


6334(6334b) 


cl 146 


I7 h 20 m 3.3 s 


-36° 9'25.9" 


0.22 


1.21 


0.33 


7.45X10 1 


7.45X10 1 


8.98x10 s 


NGC 


6334(6334b) 


cl 147 


17 h 17 m 10.9 s 


-36°26'30.8" 


0.22 


1.72 


0.47 


1.06X10 2 


1.06X10 2 


4.17x10 s 


NGC 


6334(6334c) 


cl 148 




-35°37'25.0" 


0.22 


1.82 


0.42 


1.12X10 2 


1.12X10 2 


6.57x10 s 


NGC 


6334(6334b) 


cl 149 


17 h 21 m 29.1 s 


-35°31'31.8" 


0.22 


0.67 


0.27 


4.13X10 1 


4.13X10 1 


9.08x10 s 


NGC 


6334(6334b) 


cl 151 


17' , 23 m 13.0 s 


_34°49'40.1" 


0.22 


0.99 


0.30 


6.13X10 1 


6.13X10 1 


9.83x10 s 


NGC 


6334(6334c) 


cl 152 


17 fe 17 m 41.3 s 


-36° 9T3.0" 


0.21 


1.13 


0.33 


6.97X10 1 


6.97X10 1 


8.41x10 s 


NGC 


6334(6334c) 


cl 153 


17 ,l 20 m 19.7 s 


-35°42'29.9" 


0.21 


1.17 


0.39 


7.24X10 1 


2.99X10 1 


2.18x10 s 


NGC 


6334(6334b) 


cl 155 


17 h 18 m 28.8 s 


-36°12'51.8" 


0.21 


1.29 


0.36 


7.97X10 1 


7.97X10 1 


7.40x10 s 


NGC 


6334(6334c) 


cl 156 


17 h 18 m 15.9 s 


-36°26'51.4" 


0.20 


1.14 


0.39 


7.01X10 1 


7.01X10 1 


5.12x10 s 


NGC 


6334(6334c) 


cl 159 


17 h 20 m 31.6 s 


-35°48'29.9" 


0.20 


0.46 


0.24 


2.82X10 1 


1.17X10 1 


3.65x10 s 


NGC 


&334(6334a) (6334b) 


cl 161 


17 h 19 m 58.7 s 


-36° 8'46.0" 


0.20 


0.56 


0.24 


3.48X10 1 


3.48X10 1 


1.09X10 4 


NGC 


63 34(6334b) 


cl 163 


17 ,l 22 m 14.2 s 


-35°26'49.2" 


0.20 


2.37 


0.53 


1.46X10 2 


1.46X10 2 


4.02x10 s 


NGC 


6334(6334b) 


cl 164 


17' , 16 m 13.2 s 


-36°28'25.0" 


0.20 


0.57 


0.24 


3.51X10 1 


3.51X10 1 


l.lOxlO 4 


NGC 


6334(6334c) 


cl 165 


l 7 h 18 m j 5 s 


-36°17'54.6" 


0.20 


1.21 


0.36 


7.44X10 1 


7.44X10 1 


6.91x10 s 


NGC 


6334(6334c) 


cl 167 


17 h 18 m 30.5 s 


-36°25'23.9" 


0.20 


0.82 


0.30 


5.08X10 1 


5.08X10 1 


8.16x10 s 


NGC 


63 34(6334c) 


cl 169 


17 h 20 m 48.0 s 


-35°51'49.3" 


0.19 


1.13 


0.36 


6.99X10 1 


2.26X10 1 


2.10x10 s 


NGC 


63 34(6334a) (6334b) 


cl 171 


17 h 22 m 36.2 s 


-34°57T9.4" 


0.19 


0.73 


0.30 


4.51X10 1 


4.51X10 1 


7.24x10 s 


NGC 


63 34(6334c) 


cl 175 


17 h 23 m 19.3 s 


-35°16'59.5" 


0.19 


0.86 


0.33 


5.32X10 1 


5.32X10 1 


6.42x10 s 


NGC 


6334(6334c) 


cl 176 


17 h 16 m 40.4 s 


-36°27'39.6" 


0.19 


1.66 


0.45 


1.02X10 2 


1.02X10 2 


4.86x10 s 


NGC 


6334(6334c) 


cl 179 


17 ,, 20 m 55.2 s 


-35°48'45.4" 


0.18 


0.35 


0.21 


2.16X10 1 


2.16X10 1 


l.OlxlO 4 


NGC 


6334(6334a) (6334b) 


cl 180 


17 h 17 m 37.5 s 


-36°26'24.7" 


0.18 


1.38 


0.42 


8.53X10 1 


8.53X10 1 


4.99x10 s 


NGC 


63 34(6334c) 


cl 181 


17 fc 16 m 53.0 s 


-36°27'49.0" 


0.18 


1.52 


0.45 


9.39X10 1 


9.39X10 1 


4.47x10 s 


NGC 


6334(6334c) 


cl 182 


17 fc 21 m 56.5 s 


-35°26'42.4" 


0.18 


1.35 


0.36 


8.33X10 1 


8.33X10 1 


7.73x10 s 


NGC 


6334(6334b) 


cl 184 


17 fe 21 m 5 gs 


-35°48'52.9" 


0.18 


0.64 


0.27 


3.97X10 1 


3.97X10 1 


8.73x10 s 


NGC 


6334(6334b) 


cl 185 


17^20 m 44.0 s 


-35°39'49.3" 


0.18 


0.48 


0.24 


2.95X10 1 


2.95X10 1 


9.23x10 s 


NGC 


6334(6334b) 


cl 186 


17 h 18 m 7.1 s 


-36° 7'54.8" 


0.18 


1.68 


0.45 


1.04X10 2 


1.04X10 2 


4.94x10 s 


NGC 


6334(6334c) 


cl 188 


17 h 22 "i 8.5 s 


-35°10'57.7" 


0.18 


0.49 


0.24 


3.06X10 1 


3.06X10 1 


9.58x10 s 


NGC 


6334(6334c) 


cl 189 


17 h 22 m 30.5 s 


-35°25T9.9" 


0.18 


0.46 


0.24 


2.87X10 1 


2.87X10 1 


8.98x10 s 


NGC 


6334(6334b) 


cl 190 


l 7 h l8 m 7 Q s 


-36°12'50.8" 


0.18 


0.78 


0.30 


4.83X10 1 


4.83X10 1 


7.75x10 s 


NGC 


6334(6334c) 


cl 192 


l7 h 2S m 6.1 s 


-35°13'32.9" 


0.17 


0.81 


0.33 


5.01X10 1 


5.01X10 1 


6.04x10 s 


NGC 


6334(6334c) 


cl 194 


17 h 20 m 56.4 s 


-35°32'21.1" 


0.17 


2.63 


0.53 


1.62X10 2 


1.62X10 2 


4.47x10 s 


NGC 


6334(6334b) 


cl 197 


17 h 17 m 58.4 s 


-36°12T8.4" 


0.17 


0.58 


0.27 


3.56X10 1 


3.56X10 1 


7.84x10 s 


NGC 


63 34(6334c) 


cl 198 


17 h 20 m 54.0 s 


-35°50'37.3" 


0.17 


0.59 


0.27 


3.66X10 1 


3.66X10 1 


8.06x10 s 


NGC 


63 34(6334a) (6334b) 


cl 200 


17 h 19 m 42.1 s 


-36°22'37.9" 


0.17 


1.09 


0.36 


6.74X10 1 


6.74X10 1 


6.26x10 s 


NGC 


63 34(6334c) 


cl 201 


17 h 19 m 48.2 s 


-35°47'42.0" 


0.17 


0.44 


0.24 


2.72X10 1 


2.72X10 1 


8.52x10 s 


NGC 


6334(6334b) 


cl 203 


l 7 h 20 m 2 .6 s 


-36° 0' 6.1" 


0.17 


0.96 


0.36 


5.93X10 1 


5.93X10 1 


5.51x10 s 


NGC 


6334(6334a) (6334b) 


cl 204 


17 fc 21 m 18.8 s 


-35°42T2.2" 


0.16 


0.58 


0.27 


3.58X10 1 


3.58X10 1 


7.87x10 s 


NGC 


6334(6334a) (6334b) 


cl 205 


l 7 h 18 m 5 7 s 


-36°11'30.5" 


0.16 


1.43 


0.45 


8.82X10 1 


8.82X10 1 


4.19x10 s 


NGC 


63 34(6334c) 


cl 208 


17 fc 19 m 38.1 s 


-36°23'33.7" 


0.16 


0.57 


0.27 


3.51X10 1 


3.51X10 1 


7.74x10 s 


NGC 


6334(6334c) 


cl 209 


17 fc 18 m 35.9 s 


-36°23'48.1" 


0.16 


0.88 


0.36 


5.45X10 1 


5.45X10 1 


5.06x10 s 


NGC 


6334(6334c) 


cl 210 


17 h 19 m 58.0 s 


-35°55'58.1" 


0.16 


1.34 


0.42 


8.28X10 1 


***** 


***** 


NGC 


6334(6334a) (6334b) 


cl 213 


17 h 22 m 2 A s 


-35°23'30.1" 


0.16 


1.09 


0.39 


6.74X10 1 


6.74X10 1 


4.92x10 s 


NGC 


6334(6334b) 


cl 214 


17' , 20 m 11.2 s 


-36°14T3.9" 


0.15 


0.20 


0.15 


1.24X10 1 


1.24X10 1 


1.59X10 4 


NGC 


63 34(6334c) 


cl 215 


17 h 20 m 38.8 s 


-35°43' 1.6" 


0.15 


0.55 


0.27 


3.38X10 1 


1.40X10 1 


3.08x10 s 


NGC 


Q334(6334a) (6334b) 


cl 216 


17 fe 17 m 33.2 s 


-36°13'36.5" 


0.15 


0.53 


0.27 


3.27X10 1 


3.27X10 1 


7.19x10 s 


NGC 


63 34(6334c) 


cl 217 


17 h lg m q gs 


-36°17'22.2" 


0.15 


0.87 


0.33 


5.37X10 1 


5.37X10 1 


6.48x10 s 


NGC 


6334(6334c) 


cl 218 


17 h 21 m 56.6 s 


-35°28T8.5" 


0.15 


0.92 


0.33 


5.66X10 1 


5.66X10 1 


6.83x10 s 


NGC 


6334(6334b) 


cl 219 


17' , 23 m 52.1 s 


-34°51'51.5" 


0.15 


0.49 


0.27 


3.01X10 1 


3.01X10 1 


6.63x10 s 


NGC 


6334(6334c) 


cl 220 


17 fc 21 m 17.0 s 


-35°53' 8.5" 


0.15 


0.27 


0.18 


1.68X10 1 


1.68X10 1 


1.25X10 4 


NGC 


63 34(6334b) 


cl 223 


17 fe 21 m 23.4 s 


-35°41' 0.2" 


0.15 


0.16 


0.15 


l.OlxlO 1 


l.OlxlO 1 


1.30X10 4 


NGC 


63 34(6334a) (6334b) 


cl 226 


17' , 22 m 38.5 s 


-35°29T9.3" 


0.15 


0.44 


0.24 


2.70X10 1 


2.70X10 1 


8.46x10 s 


NGC 


6334(6334b) 


cl 227 


17 ,, 20 m 46.8 s 


-35°59'25.4" 


0.15 


0.57 


0.27 


3.50X10 1 


3.50X10 1 


7.70x10 s 


NGC 


6334(6334b) 


cl 229 


17 h 20 m 52.7 s 


-35°59'49.2" 


0.14 


0.49 


0.27 


3.03X10 1 


3.03 xlO 1 


6.68x10 s 


NGC 


6334(6334b) 


cl 231 


17 h 23 m 11.4 s 


-35°14'52.4" 


0.14 


0.83 


0.33 


5.12X10 1 


5.12X10 1 


6.17x10 s 


NGC 


6334(6334c) 


cl 232 


17 fc 21 m 47.7 s 


-35°39'47.2" 


0.14 


0.32 


0.21 


1.95X10 1 


1.95X10 1 


9.10x10 s 


NGC 


6334(6334b) 


cl 233 


17 ,, 23 m 13.3 s 


-35°15'32.0" 


0.14 


0.76 


0.33 


4.69X10 1 


4.69X10 1 


5.66x10 s 


NGC 


6334(6334c) 


cl 237 


17 h 19 m 16.5 s 


-36° 2' 5.6" 


0.14 


0.80 


0.33 


4.91X10 1 


4.91X10 1 


5.92x10 s 


NGC 


6334(6334b) 


cl 240 


17 h 20 m 57.9 s 


-35°49'57.0" 


0.14 


0.46 


0.24 


2.83X10 1 


2.83X10 1 


8.87x10 s 


NGC 


6334(6334a) (6334b) 
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cl 241 


\7 h 22 m 7.6 s 


-35°25'29.6" 


0. 


14 


0.41 


0.24 


2.52X10 1 


2.52X10 1 


7.90x10 s 


NGC 


6334(6334b) 


cl 242 


17 fc 18 m 17.6 s 


-36° 9'47.2" 


0. 


14 


0.34 


0.21 


2.12X10 1 


2.12X10 1 


9.94x10 s 


NGC 


6334(6334c) 


cl 243 


17 h 19 m 56.7 s 


-35°53'42.0" 


0. 


14 


0.41 


0.24 


2.53X10 1 


2.53X10 1 


7.94X10 3 


NGC 


6334(6334a) (6334b) 


cl 245 


l7 h 20 m 197 s 


-35°48'54.0" 


0. 


14 


0.49 


0.27 


3.03X10 1 


1.25X10 1 


2.75X10 3 


NGC 


6334(6334a) (6334b) 


cl 246 


17 fe 19 m 51.4 s 


-36° 2'30.1" 


0. 


14 


0.62 


0.30 


3.80X10 1 


3.80X10 1 


6.11X10 3 


NGC 


63 34(6334b) 


cl 247 


17 fe 17 m 32.8 s 


-36°26' 8.5" 


0. 


13 


0.66 


0.30 


4.05X10 1 


4.05X10 1 


6.49X10 3 


NGC 


633A(6334c) 


cl 249 


17' l 22 m 19.3 s 


-35°21'20.9" 


0. 


13 


1.42 


0.45 


8.76X10 1 


8.76X10 1 


4.17x10 s 


NGC 


6334(6334b) 


cl 251 


17 ,i 23 m 36.5 s 


-34°51'45.4" 


0. 


13 


0.32 


0.21 


2.00X10 1 


2.00X10 1 


9.36X10 3 


NGC 


6334(6334c) 


cl 252 


17 h 19 m 40.2 s 


-36° 3' 9.7" 


0. 


13 


0.39 


0.24 


2.38X10 1 


2.38X10 1 


7.45x10 s 


NGC 


6334(6334b) 


cl 253 


17' i 20 m 50.1 s 


-35°56'37.3" 


0. 


13 


0.43 


0.24 


2.64X10 1 


2.64X10 1 


8.27X10 3 


NGC 


6334(6334b) 


cl 256 


17' , 20 m 53.4 s 


-36° 1'57.4" 


0. 


13 


0.45 


0.27 


2.80X10 1 


2.80X10 1 


6.16X10 3 


NGC 


6334(6334c) 


cl 259 


17 h 21 m 6.1 s 


-35°23'33.0" 


0. 


13 


0.09 


0.12 


5.56 


5.56 


1.39X10 4 


NGC 


6334(6334c) 


cl 260 


l 7 h l8 m 9A s 


-36°20T9.0" 


0. 


13 


0.72 


0.33 


4.43X10 1 


4.43X10 1 


5.34X10 3 


NGC 


633A(6334c) 


cl 261 


17 fe 21 m 17.7 s 


-35°51'56.5" 


0. 


13 


0.38 


0.24 


2.33X10 1 


2.33X10 1 


7.30X10 3 


NGC 


6334(6334b) 


cl 262 


17' , 20 m 35.6 s 


-36° 3' 1.8" 


0. 


13 


0.19 


0.18 


1.19X10 1 


1.19X10 1 


8.81X10 3 


NGC 


6334(6334b) 


cl 264 


17 fe 19 m 33.5 s 


-36°24'37.8" 


0. 


12 


0.21 


0.18 


1.28X10 1 


1.28X10 1 


9.50X10 3 


NGC 


6334(6334c) 


cl 266 


I7 h 22 m 2.1 s 


-35° 13' 5.9" 


0. 


12 


0.34 


0.24 


2.13X10 1 


2.13X10 1 


6.68X10 3 


NGC 


6334(6334c) 


cl 268 


17' , 21 m 12.7 s 


-35°30'36.7" 


0. 


12 


0.78 


0.36 


4.83X10 1 


4.83X10 1 


4.49 xlO 3 


NGC 


6334(6334b) 


cl 271 


17 h 22 m 25.1 s 


-35°16'56.3" 


0. 


12 


0.79 


0.36 


4.90X10 1 


4.90X10 1 


4.55X10 3 


NGC 


633A(6334c) 


cl 273 


17 h 21 m 25.8 s 


-35°31'32.2" 


0. 


12 


0.59 


0.30 


3.66X10 1 


3.66X10 1 


5.88X10 3 


NGC 


633A(6334b) 


cl 276 


17 h 16 m 57.0 s 


-36°25'25.3" 


0. 


12 


0.38 


0.24 


2.37X10 1 


2.37X10 1 


7.41x10 s 


NGC 


6334(6334c) 


cl 279 


17 h 2Q m 5.3 s 


-36°27T0.1" 


0. 


12 


0.22 


0.18 


1.36X10 1 


1.36X10 1 


l.OlxlO 4 


NGC 


6334(6334c) 


cl 282 


17' , 20 m 12.6 s 


-36° 7' 1.9" 


0. 


11 


0.17 


0.15 


1.06X10 1 


1.06X10 1 


1.36X10 4 


NGC 


6334(6334b) 


cl 284 


17 ,1 22 m 22.1 s 


-35°27' 4.7" 


0. 


11 


0.73 


0.36 


4.51X10 1 


4.51X10 1 


4.19x10 s 


NGC 


633A(6334b) 


cl 286 


17 h 19 m 42.9 s 


-35°59'49.9" 


0. 


11 


0.44 


0.27 


2.74X10 1 


2.74X10 1 


6.02x10 s 


NGC 


6334(6334b) 


cl 287 


17' i 23 m 37.8 s 


-34°52'57.4" 


0. 


11 


0.27 


0.21 


1.67X10 1 


1.67X10 1 


7.80x10 s 


NGC 


6334(6334c) 


cl 290 


l 7 h 2Q m 7 3 s 


-36° 7'34.0" 


0. 


11 


0.18 


0.18 


1.14X10 1 


1.14X10 1 


8.44x10 s 


NGC 


6334(6334b) 


cl 292 


17' , 23 m 48.8 s 


-34°52T6.0" 


0. 


11 


0.42 


0.24 


2.58X10 1 


2.58X10 1 


8.09x10 s 


NGC 


6334(6334c) 


cl 293 


17' , 22 m 13.1 s 


-35°35' 5.3" 


0. 


11 


0.14 


0.15 


8.40 


8.40 


1.08X10 4 


NGC 


633A(6334b) 


cl 294 


l 7 h 23 m 9 7 s 


-34°49'32.2" 


0. 


11 


0.51 


0.27 


3.14X10 1 


3.14X10 1 


6.92x10 s 


NGC 


6334(6334c) 


cl 295 


17 ,i 22 m 17.3 s 


-35°19'53.0" 


0. 


11 


1.43 


0.47 


8.86X10 1 


8.86X10 1 


3.47x10 s 


NGC 


6334(6334b) 


cl 300 


17 h 17 m 44.2 s 


-36°22'33.2" 


0. 


11 


0.19 


0.18 


1.20X10 1 


1.20X10 1 


8.90x10 s 


NGC 


6334(6334c) 


cl 303 


17 fc 17 m 48.2 s 


-36°22T7.4" 


0. 


11 


0.24 


0.21 


1.49X10 1 


1.49X10 1 


6.99x10 s 


NGC 


6334(6334c) 


cl 307 


17 h 23 m 3.4 s 


-35°12'21.2" 


0. 


11 


0.14 


0.15 


8.71 


8.71 


1.12X10 4 


NGC 


6334(6334c) 


cl 308 


17 h 23 m 5.8 s 


-35° 5'33.0" 


0. 


11 


0.36 


0.24 


2.24X10 1 


2.24X10 1 


7.01x10 s 


NGC 


6334(6334c) 


cl 309 


17 h 18 m 16.1 s 


-36°18'35.3" 


0. 


11 


0.88 


0.39 


5.42X10 1 


5.42X10 1 


3.96x10 s 


NGC 


633A(6334c) 


cl 311 


17' i 22 m 35.2 s 


-35°28'55.6" 


0. 


11 


0.21 


0.18 


1.28X10 1 


1.28X10 1 


9.50x10 s 


NGC 


6334(6334b) 


cl 318 


17 fe 17 m 15.8 s 


-35°59'51.0" 


0. 


10 


0.06 


0.09 


3.95 


3.95 


2.35X10 4 


NGC 


6334(6334c) 


cl 319 


17 fe 16 m 58.0 s 


-36°17'49.6" 


0. 


10 


0.87 


0.36 


5.36X10 1 


5.36X10 1 


4.98x10 s 


NGC 


6334(6334c) 


cl 326 


17 fc 21 m 34.5 s 


-35°37'39.7" 


0. 


10 


0.12 


0.15 


7.16 


7.16 


9.20x10 s 


NGC 


633A(6334b) 


cl 329 


17 20 24.2 


-35 2y 25.8 


0. 


10 


0.17 


0.18 


1.06X10 1 


l.Obx 10 


7.89x10 s 


NGC 


633A(6334c) 


cl 336 


17 h 19 m 42.8 s 


-36° 3' 1.8" 


0. 


10 


0.28 


0.21 


1.73X10 1 


1.73x 10 1 


8.09x10 s 


NGC 


6334(6334b) 


cl 337 


17 h 19 m 42.8 s 


-36°23'25.8" 


0. 


10 


0.29 


0.21 


1.80X10 1 


1.80X10 1 


8.41x10 s 


NGC 


6334(6334c) 


cl 340 


l 7 h 20 m g gs 


-36° 6'38.2" 


0. 


10 


0.13 


0.15 


8.21 


8.21 


1.05X10 4 


NGC 


6334(6334b) 


cl 342 


17 fc 21 m 25.9 s 


-35°32'36.2" 


0. 


10 


0.17 


0.15 


1.05X10 1 


1.05X10 1 


1.35X10 4 


NGC 


6334(6334b) 


cl 343 


l 7 h 22 m 2QA s 


-35°35T2.8" 


0. 


10 


0.20 


0.18 


1.27X10 1 


1.27X10 1 


9.41x10 s 


NGC 


633A(6334b) 


cl 344 


17 ,l 20 m 25.0 s 


-35°41T0.0" 


0. 


10 


0.21 


0.18 


1.28X10 1 


1.28X10 1 


9.54x10 s 


NGC 


633A(6334b) 


Note. 


— The clump 


numbers highlighted 


in 


the first 


column 


arc those which are 


considered 


to be significantly warmer. Clump 210 has 


undetermined mass after 


subtracting the frce- 


-free emission at 1.2 


mm 













a Afi.2mm is the computed clump mass from all the emssion at 1.2 mm and assuming a temperature of 17 K, while M considers the two- 
tcmpcraturc clump ensemble and the correction by free-free emission. 

k The average number density was computed using i? c ff and M and assuming spherical clumps with a mean molecular weight of (X — 2.3. 



